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THE DETERMINATION OF THE QUANTITY OF RADIOACTIVITY 
CONTAINED IN A GASEOUS EFFLUENT AT THE NEVADA TEST SITE 
USING GAMMA DOSE-RATE MEASUREMENTS 


I. Introduction 


Effluent documentation may be simply stated as the determination of the 
radioactive materials contained in an effluent release. The quality of the effluent 
is usually determined by detailed laboratory analyses, whereas the quantitative 
determinations are the results of relatively simple calculations made at the 
time of release. This paper will concern itself only with the quantitative deter¬ 
minations presently used to arrive at effluent release estimates. 


II. Methods 

Several different geometric configurations of an effluent cloud are possible 
with respect to a radiation detector. These may be the result of different escape 
velocities, meteorological conditions, wind speeds and directions, relative dis¬ 
tances, relative dimensions, and angles, etc. 

Uncertainty as to possible time of appearance of an effluent and lack of 
knowledge of wind speed and direction which may possibly exist at the time of 
appearance make it judicious that preparations include a concentric array of 
radiation detectors equally spaced around the potential effluent release point(s). 

The distances and spacings of the radiation detectors are arbitrary. A 
typical array is shown below in which eight detectors are placed on a 1000-foot 
diameter circle and eight detectors are placed on a 2000-foot diameter circle 
around the central potential release point. 
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A. Line Sources 

1. Vertical Line Source 

If the gaseous effluent assumes a geometric configuration with respect to 
the array of selectors, similar to that shown in the diagram below, the applicable 
equations used to arrive at the quantity of effluent released may be derived as 
follows: 
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The differential photon flux at point P from differential source dx, is 


S Be "^^dx 

= - - 

2 

An r 


(Reference 1) 


where 


_ / photons 

(j) = flux at point of measurement P I ~ 


mm sec 


S = specific line source strength( - ■ 

X \cmseCi 


(1) 


B(r) = build-up factor 

r = slant distance from increment of x, dx, to point of 
measurement P (cm) 

X = distance from ground to source increment, dx (cm) 

= total attenuation coefficient (cm ^) 
h = source height (cm) 

If we assume a uniform line source strength, S is then constant, and 


mhotons ^ 
l sec j 


= S 


X 


photons ^ 
cm sec ; 


h(cm) 


( 2 ) 
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For gross fission products. 


S 


^ photons ] 
^ sec / 


Q (curies) • 



(Reference 2) 


For gross fission products at early times up to approximately H + 30 min¬ 
utes, the average gamma photon energy used in these calculations is ~ 0. 92 5 

(References 2, 4). Corrections may be easily made for any other energy 
photon ' 

desired. For this particular energy, 1 roentgen per hour is approximately 

5. 7 X 10^ (Reference 3). However, the effluent release may be highly 

cm sec ^ . 

fractionated and obviously the composition of any release changes with time; 

therefore, the appropriate and conversion factor must be selected in each 

particular case. 

Combining Equations 2 and 3, solving for S , and substituting Equation 1, 
we have: 


3.7xl0^“( -Q(curies)dx 

_ \curie sec/ _ \ dis / __ 

2 2 ^ 
477-h r (cm ) 

Since we are assuming a uniform concentration, e.g., a constant linear source 
strength: 


d(/)( 


ihotonsN 


\cm sec 




Q _ dQ . - Q 
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Converting the flux to a gamma dose rate: 


dcj) 


roentge n s \ 
\ hour / 


photons' 


mm sec/ 


( 7 ) 


5. 7 X 10^ 


' photons hour 
roentgens cm^ sec 


jp /roentgens j 
\ hour / 


3.7 X 10 


10 / dis 


> curie sec 


j • 1.19 (curies) 


2 , 2, ^ , ^5 / photons hour 

4r7 r (cm ) ♦ 5.7 x 10 f— - - - 

\ roentgens cm sec/ 


( 8 ) 


R 


roentgens \ =6.15x10^ / ^ 


dQ 


hour 


) 




r cm 


curies 


( 9 ) 


This general equation appears again and again in effluent release calcu¬ 
lations. To estimate Q, the total release, from this integral equation, one 
must delineate the dimensions of the source cloud, and make assumptions about 
B, and /i, and how dQ varies with r. With these assumptions, it is then possible 

to relate Q to R. 

From the diagram on page 6 

2 2 2 

r = a sec 0 r = a sec 0 

2 

X = a tan 0 dx = a sec 6 dO 
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Substituting for dx and solving for dQ in Equation 4 and substituting for r and 
dQ in Equation 8, we have: 



where we have assumed B(r) to be slowly varying, i.e., essentially constant- 


/roentgens j 
\ hour / 


6. 15 X 10 B Q (curies) 
h(cm) a (cm) 



max 


e-,-asece 


( R/hr cm 
curies 


( 11 ) 


Thus 


Q(curies) = 


R( ro^nt£ens\ _ (cm) • a (cm) 
hour / 


curies 


6.15x IC^ • B •/ e 
6=0 


rS 

I max 

7 ^ 


-/iasect?^^ ^R/hrcm^ 


( 12 ) 


where R is the dose rate reading of the gamma radiation detector at the time 
of measurement. 
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2. Symmetric Line Source 

If the line source configuration is such that it is similar to the following 
diagram. 
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the expression used for the calculation is: 


(/> ^ 


B S 


r 


0^ or 


- aa. sec ^ _ 
e ^ dd 


277a 


(Reference 1) 


(15) 


where 


X 


= specific source strength ^-— - " —), assumed constant 
^ Vcm sec / 


X = source length (cm) 


B = build-up factor 


^ = total attenuation coefficient (cm ^) 


a = perpendicular distance from source to point of measurement P (cm) 

, photons \ 

<j) = flux at point of measurement P/- - -j 

\cm sec j 


3. Uniform Release Rate into a Line Source 

If a uniform rate of release and a uniform concentration are assumed, 
the source length X can be thought of as the gaseous effluent being evenly dis¬ 
tributed over a distance directly proportional to the wind speed and the time 
duration of effluent release; e.g., an increment of x, dx, contributes a given 
flux at the point of measurement P for a given time as it traverses the distance 
from the point of release to some distance beyond which its contribution to the 
total flux at point P is insignificant. This may be represented simply by: 


12 


X (cm) 



T (sec) 


(16) 


Q(ciiries) 


S 

X 


[curies ’ 


cm 


X (cm) 


(17) 


<l> 


.*.Q( curies) 


photons 


mm secy 


_ cm 

Ana (cm) • U gg^ • T (sec) 


B 


rO 

I m: 

. 

min 


max _ sec0 

e ^ dd 


(18) 


where 


U = average wind speed 
T = time duration of effluent release (sec) 



By the same reasoning and the same conditions as previously used. 


5.7 X 10^ _ photons hour ] ^ ^ (roent gen s ^ 4^a(cm). ul.£mj. T(sec) 


Q(curies) = 


^roentgens cm sec/ 


hour 


iSecy 


3.7 X 10 


10 


dis \ -j^ iQ photon^ ^ f^ax^-^asecd 


murie seel 


dis 


de 


(19) 


6 . 

min 


Q(curies) = 


1.63xl0“tRp|B^j.a(cm). T(sec) 


d . 

min 


-ga sec0 





curies 


cm 

cm -- 

sec 


sec} 


( 20 ) 
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Converting from metric to English units. Equation 20 becomes 


Q(curies) = 


1. 63 X lo'S Roentgens' ! . . jj- 

V hour / 


— / cm' 
seel 


T (sec) 


B 


r^a 

^in 


-iiasec0,„ 
e ^ dd 


curies 


tR/hr cm"^ 


4.47 X 10^ cm/sec ^ 3. 05 x 10^ cm ^ 6 x 10^ sec 
1 mile /hour ft min 


( 21 ) 


Q(curies) = 


13.2 R / roentgens ) . ^ fftl • Utiles • t fmini 
' hour / hour 


B 


/ 


6 


-as. seed 
e ^ dd 


mm 


Again if 


curies 


iR/hr • ft • • mini 


hour 


( 22 ) 


then 


d. 


min' 


max' 



min 


-a a. seed 
e ^ 


dd 



0 . or 

min' 


0 1 

-u a. seed.a 
e dd 
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and the equation becomes 


Q(curies 


6.6R /i:2entHerm\ 

_\ hour / 


a (ft) • U 


/ miles\ 
\hour / 


T (min) 


B 


I 


(23) 


^min ^max 


e -pLcL seed 


dd 


The factor a is determined by use of maps, drawn to scale, and the 
average value of the wind direction over the time period for which the calcula¬ 
tions are made. The following diagram shows the relation of a detector to a 
line source for a given wind direction. 



By use of ruler and protractor, the approximate position of the line source can 
be determined. The positioning of this line is arbitrary; however, the final 
values of Q determined by this method are usually an average of two or more 
detectors each having its own particular values for the parameters in the equa¬ 
tion at the time the measurements are made. 
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Whenever possible optical measurements are obtained of the effluent 
cloud, thus giving more reliable values of dimensions and relative locations. 
However, the effluent is often not visible, nor does the visible cloud always 
necessarily correspond to the areas or locations of the radioactive effluent 
c one ent r at ion s. 

The a measured above may not necessarily be the effective source to 
detector distance. Diffusion of the effluent during the time of traversal of the 
distance from the point of release to the array of detectors may result in an 
effective source center of activity, assuming uniform concentration, and thus 
an effective source distance. 



The extent of diffusion, of course, depends on the different meterological param¬ 
eters. Figure 1 represents an attempt to simplify these determinations such 
that the effective source distance may be arrived at more rapidly. 

The resToltant figures are based entirely on theoretical considerations 
using those values for c and n for the Nevada Test Site shown on Figure 1. 

The function e must be evaluated for the different 

^min 

values of a and 6. These data are plotted in Figure 1-9, page 385 of Reference 1. 

Applying this function to a gaseous fission product effluent, using the average 

- 5 -1 

photon energy as 0. 925 mev and a total attenuation coefficient of 8.19 x 10 cm 
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Figure 1. 


Nomograph for Determination of Half-Height of Cloud 





_ _ 3 _ _ 

(Reference 2), with the values of 1 meter < a < 10 meters and 20° < 0 -v 90°, 
results in the data given in Table I. For convenience in its use for the solution 


of the preceding equations, ^ —-- — - 

/ “ LLSL SQC0 

e ^ d^^ versus; 

for 6 = 90°. These data are plotted in Figure 3. 

B. Point Source 

For an isolated quantity of radioactivity in a geometrical configuration 
such that the dimension of the source is small compared to the distance from 
the source to the point of measurement, a point source equation may be used 
for the calculation. 


is plotted versus 6 for the selected 


(Cloud) 


B S e 
o 


(Reference 1) 


where 


<j) = flux at point of measurement P 


photons 


cm sec 


B = build-up factor 

S = source strength 
o 

M = total attenuation coefficient (cm ") 

r = distance from source to point of measurement P (cm) 
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TABLE I 


Values of F(0 


.b) =f 

Or^ 


'max -/xasec^,., , 

e b = \ia., 


^min 

for selected values of a and $ 




when ii = 

8.19 X 10 m 



a 

_ 

F(0.b) 

[6 = 20°) 

no. b) 

(6 = 40°) 

F(0,b) 

(d = 60°) 

F(0,b) 

(0 = 90°) 

10 m 

8. 19 X lo'^ 

3.2 X lO”^ 

6. 4 X lO"^ 

9.5 X lO"^ 

1. 3 

25 m 

2.05 X lO"^ 

2. 8 X 10'^ 

5. 6 X lO"^ 

8.0 X 10'^ 

1. 0 

50 m 

4.1 X 10“^ 

2. 3 X lO"^ 

4. 4 X lO"^ 

6.2 X lO"^ 

7.2 X lO"^ 

100 m 

8.19 X lO'^ 

1.4 X lO"^ 

2. 8 X lO'^ 

3.8 X 10“^ 

4.1 X lO"^ 

200 m 

1. 64 

6.4 X lO"^ 

1. 15 X lO"^ 

1.4 X lO'^ 

1. 4 X lO'^ 

300 m 

2. 46 

2.8 X lO"^ 

4. 8 X lO"^ 

5.6 X lO'^ 

5. 6 X 10"^ 

400 m 

3.28 

1.23 X 10‘^ 

2. 0 X lO"^ 

2.2 X lO"^ 

2.2 X lO"^ 

500 m 

4. 10 

5. 3 X lO"^ 

8.4 X 10‘^ 

9.2 X lO"^ 

9.2 X lO"^ 

600 m 

4. 91 

2. 3 X lO'^ 

3. 5 X 10‘^ 

3.8 X lO"^ 

3. 8 X lO"^ 

800 m 

6. 55 

4. 3 X lO"^ 

6.2 X lo"^ 

6.4 X lo"^ 

6. 4 X 10"“^ 

1000 m 

8, 19 

8.2 X 10‘^ 

1.3 X 10“^ 

1.3 X lO"^ 

1. 3 X lO"^ 
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Figure 2. 


f -fia sec6,^ 


Versus 0 for Selected Values of a 








TABLE II 

/ -ua seed ,^ , 

e d0, b = ua, 

for selected values of a when 

-3 -1 

6 = 90° and = 8. 19 x 10 m 


a (feet) 

F(<9, b) 

a (feet) 

F(0. b) 

a (feet) 

F(0, b) 

1 

1. 5 

1100 

4 X lO"^ 

2500 

9 X lO""^ 

25 

1.35 

1200 

3. 1 X lO"^ 

2600 

6. 8 X lo”'^ 

50 

1, 18 

1300 

2. 3 X lO"^ 

2700 

5.2 X 10“^ 

75 

1.08 

1400 

1.75 X lO'^ 

2800 

-4 

4x10 

100 

9. 4 X lO'^ 

1500 

1.35 X lO'^ 

2900 

3.1 X lO"^ 

200 

6.2 X 10“^ 

1600 

1.2 X lO'^ 

3000 

2.4 X 10““^ 

300 

4. 5 X 10 ^ 

1700 

7.8 X lO'^ 

3100 

1. 8 X lo"^ 

400 

3. 3 X lO'^ 

1800 

6 X lO"^ 

3200 

1.4 X 10‘‘ 

500 

2. 35 X lO"^ 

1900 

4.5 X 10~^ 

3300 

1.05 X lO" 


-1 





600 

1.75 X 10 

2000 

3. 4 X 10 



700 

1.25 X lO"^ 

2100 

2. 6 X lO'^ 



800 

9.5 X lO'^ 

2200 

2 X lO'^ 



900 

7 X lO'^ 

2300 

1.5 X lO'^ 



1000 

5. 4 X lO"^ 

2400 

1,15 X lO"^ 
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conditions as previously used. 


By the same reasoning and the same 


Q(curies) = 



(25) 


Q(curies) 


1.63 xio"^ . R 1 

roentgens \ 

, hour / 


2 2\ 
9.29x10 cm V 

B . 

e-^r 

\ 2 
\R/hr cm 

/ 


(26) 


Q{curies) = 


1.51xl0'^R(i:2™‘S£Ss). r^(tth 
\ hour /_ . 


B • e 


-fir 


curies ' 

ilVhrft2j 


(27) 


The function - or e^^ is plotted versus r in Figure 4 to facilitate the deter- 

Q-fir 

mination of Q. 

The responses of selected detectors are continuously recorded such that 
average radiation dose rates over the time of effluent release may be obtained. 

The wind speed and direction are also continuously recorded. In addition 
to this information, surveyors are also available to provide physical measure¬ 
ments of visible effluents. Cloud dimensions of height and width versus time 
are obtained and relayed via telephone,and this information often provides the 
basis for selection of the model used in arriving at the estimate of activity con¬ 
tained in the effluent. 

Again, the values used in these equations apply only to gross fission 
products for a particular period of time. For fractional releases, the values 
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used depend on the constituents of the effluent and these values will change with 
time. 


III. Activities of Certain Effluent Constituents 

The nuclides of frequent interest are the isotopes of iodine, xenon, and 
krypton. The activities of particular isotopes of these nuclides contribute pre¬ 
dominantly to certain types of effluent release. Therefore, it is of importance 
to know the maximum possible activities which might possibly be present at a 
given time. 

The activity curves for each of the nuclides shown in Figure 5 were derived 
from the data given by Bolles and Ballou (Reference 5). The listed activities as 
a function of time were converted to curies per kiloton as follows: 


(Ref. 5) (Ref. 6) 

activity (dis/sec) . i, 45 x 10 ^^ /fissions\ 


^ /curies \ 


10 fissions 


i kiloton 


3_,^10l0/dis/pc\ 
\ curie / 


(28) 


, /curies\ o 'Xo 1 ^ 

(kuStoi) ' 


8 ^ activity (dis/sec) 


10 fissions 


(29) 


These data represent only those activities resulting from thermal fission 
of U-235. Frequently, comparisons of yields of these nuclides from other types 
of fission processes must be considered. 
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Selected Isotopes of the Xenon, Krytron, 
V and Iodine Nuclides Plotted Versus Time 
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13 5 m 
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Per 

Kiloton 

of U-235 

Thermal 

Fission 


,133 




133 m 


TIME 
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Curves similar to those in Figure 5 may be prepared for the various 
fission processes. Presently, however, comparisons of total chain and inde¬ 
pendent fission yields may be made for the listed nuclides by use of Table III 
(Reference 7). This table gives the total number of atoms of each nuclide 

4 

resulting from 10 fissions of each process. 


IV. Ventilation Exhaust Calculations 

When an effluent is discharged in a controlled manner, particulate and 
charcoal filters are employed. The resultant constituents of the gaseous dis¬ 
charge are the isotopes of xenon and krypton. The quantity of these gases which 
is released over a period of time may be determined by use of the following 
equations. 

The radiation dose-rate reading of the detector is proportional to the 
energy fluence rate resulting from a given volume of radioactive gases emitting 
the radiations. 



(Ji£^ . 3.7 X nn^° / dis ' . to-6/ourle\ 
\cm^/ \curie sec/ \ j 


1 1 p / photon\ 

. 1 mev 

\ dis / 

photon 


"b 7 mev hour ^(Ref. 3) 

5. 5 X 10 -- 2 -1 

I photons cm sec/ 


.-Mr 


dV (cm ) 


2 2 
4n r (cm ) 


(30) 
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TABLE III 


Total Number of Atoms of Selected Nuclides Resulting from 10 Fissions 
of Several Different Fission Processes (Reference 7) 



U235 

TO 239 

Pu 

U235 

239 

Pu 

U235 

T3 239 
Pu 


(thermal) 

(thermal) 

(14 mev) 

(14 mev) 

(Fast) 

(Fast) 

Xe^2^ 


1.3036 

3.1899 

8.5891 


1.0964 

Xe'3' 

93.0372 

135.3368 

132.4311 

165.6841 

54.8015 

125.8819 


275.1100 

288.5859 

210.7298 

210.4086 

300.9506 

280.1888 

T. 85 

Kr 

1.1104 

1.3335 

9.8598 

6.1293 

1.0371 

1. 7258 

87 

Kr 

32.8131 

18.8078 

67.7622 

36.0846 

32.9419 

29.1356 

CO 

CO 

86. 9348 

43.1884 

127.5979 

68.3066 

83.6765 

44.0587 

ll31 

0.9454 

4.6273 

12.6631 

17.8898 

1.1794 

5.6139 

jl32 

16.7067 

29.8537 

41.5660 

50.4168 

19.4194 

29.8462 

jl33 

64.4258 

91. 6509 

93.4751 

107.7455 

54.2931 

76.8489 

jl35 

285.1140 

329.3860 

260.9957 

263.2649 

236.8970 

304.5531 
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/roentgens \ 
\ hour / 


= 6. 4 X 10 



dV , , 

— 2 - (cm) 



(31) 


The factor of e"^^ is considered to be a relatively insignificant correction 
under these particular conditions since other uncontrollable factors contribute 
greater errors in the calculations. The value of = 1 mey was used purpose¬ 
fully in this expression. As the relative composition of the effluent changes with 
time, the effective energy changes and thus it is a simple matter to make the 
energy corrections when necessary. 

A radiation detector is placed adjacent to a section of exhaust pipe as 
shown below: 


dx 



An approximate solution to the function J in terms of this diagram 

r2 




is as follows: 


V = ttD X 


(32) 


r^ = (D + d)^ + 


dV = ttB dx 


(33) 

(34) 
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=/— 


(D + d)^ + 


.- = C 

d)^ + X^ 


(D + d)^ + x^ 


jy = I ]__ arc tan f . . . | + j -~rT TTTTr 

^ \D+d D+d/ \D+d TTTT 


2 / ^ 

dV = nP j arc tan ^— + arc tan - 

r^ D+d\ D+d D+d 


Equation 31 may be rearranged and converted to English units. 


{ 5 )- 


roentgens 
^ hour i 


6.4xl0'^*/^(cm) \R/hrcm3 10 mc 


4 3 

ftc cm curie 2> 8 x 10 cm 

/ Vi-yy 3 1 ■P+^ 


( _^_ \ 

\3. 05 X 10^ cm / 


.ifcuries\ 
<—) = 


1.45x 10-^R (roentge ns 
\ hour 


/f<«) 


curies 
^R/hr ft*^ 
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When effluent is being discharged, the activity exhausted is simply: 


1.45 X 


Q(curies) 


1 / roentgen s) _w(jii).T (min) 

\ hour /_\min /_ 


/ 


^ (ft) 

r 


curies 


vR/hr ft"" 


(40) 
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